Competing Orders in Coupled Luttinger Liquids 
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We consider the problem of two coupled Luttinger liquids both at half filling and at low doping 
levels, to investigate the problem of competing orders in quasi-one-dimensional strongly correlated 
systems. We use bosonization and renormalization group equations to investigate the phase di- 
agrams, to determine the allowed phases and to establish approximate boundaries among them. 
Because of the chiral translation and reflection symmetry in the charge mode away from half filling, 
orders of charge density wave (CDW) and spin-Peierls (SP), diagonal current (DC) and d-density 
wave (DDW) form two doublets and thus can be at most quasi-long range ordered. At half-filling, 
Umklapp terms break this symmetry down to a discrete group and thus Ising-type ordered phases 
appear as a result of spontaneous breaking of the residual symmetries. Quantum disordered Haldane 
phases are also found, with finite amplitudes of pairing orders and triplet counterparts of CDW, SP, 
DC and DDW. Relations with recent numerical results and implications to similar problems in two 
dimensions are discussed. 

PACS numbers: PACS: 71.10.Fd, 71.10.Hf, 71.30.+h, 74.20. Mn 



I. INTRODUCTION 

The problem of the nature of the phase diagram of 
the cuprate superconductors remains at the center of 
research in the physics of strongly correlated electron 
systems. Recent work has focused on the possible com- 
peting orders responsible for the known features of the 
phase diagram as well as to the unusual physical prop- 
erties of the pseudo-gap regime. A number of candidate 
competing orders have been considered, including anti- 
ferromagnetism (AF), d-wave pairing (DSC), incommen- 
surate charge ordered states and other liquid crystal-like 
phases, and d-density wave states (DDW) (also known as 
staggered flux states (SF) or orbital antiferromagnetism 
(OAF)), among others. 

50(5) theory Q focuses on the competition between 
antiferromagnetism and d-wave superconductivity. In 
this theory the natural SU(2) x U(l) symmetry of the 
spin and charge degrees of freedom is regarded as the 
result of an explicit symmetry breaking of a larger sym- 
metry, characterized by a global SO (5) group. In this 
picture, this larger symmetry is not apparent except close 
to a quantum critical point whose quantum fluctuations 
suppress both antiferromagnetism and d-wave supercon- 
ductivity, thus leading to a pseudo-gap regime controlled 
by this fixed point. 

In contrast, in the stripe mechanism |2(, the ground 
state of the doped Mott insulator is an inhomogeneous 
charge ordered state resembling a liquid crystal phase , 
which breaks both rotational invariance and (partially) 
translation invariance, i. e. it is a quantum smectic. In 
this picture the pseudo-gap is the spin gap which devel- 
ops in these quasi-one-dimensional states, and it is not a 
signature of some sort of long range order. In this mech- 
anism, macroscopic phase coherence and d-wave super- 
conductivity result from inter-stripe Josephson couplings 



In the d-density wave state, and similarly in the physi- 
cally equivalent staggered flux and orbital antiferromag- 
netic states, there is a hidden order which has the same 
d x 2_ y 2 symmetry as a d-wave superconductor. In this 
phase the ground state has an ordered pattern of stag- 
gered orbital currents, and this is the order which com- 
petes with d-wave superconductivity H H 0H H E3 • 

However, in spite of a continued effort during the past 
decade or so, and largely due to the lack of systematic 
non-perturbative methods in two dimensions, it has been 
quite difficult to establish the phase diagram of reason- 
able two dimensional strongly correlated systems based 
on the Hubbard model. Much of the work done is based 
on mean field-type approximations which favor one type 
of order over others or privileges the competition among a 
particular pair of order parameters. While it is quite pos- 
sible that these studies reveal different aspects of possi- 
ble phase diagrams of some generic, possibly short range 
models, it is not possible at present to determine reliably 
the phase diagram of many of these models except some- 
times at extreme regimes of some parameter. Thus dif- 
ferent approaches, including large- N methods (and their 
relatives), have been used to construct spin-liquid states 
[H El El El El 13 113 Hartree-Fock, large-d and 
large N methods have been used to study phase sep- 
aration and striped states 0, 0, |2l [2]J. Similarly 
Hartree-Fock methods have also been used to study the 
competition between superconductivity and DDW order 
|22j| . There is also an extensive literature on numerical 
simulations which work either at moderate to high tem- 
peratures (as in Quantum Monte Carlo simulations due 
to the fermion sign problem) or exact diagonalizations 
of systems which are usually too small to resolve these 
issues. 

It is largely for these reasons, as well as for the need 



of non-perturbative results, that some of these questions 
have been considered in the framework of quasi-one di- 
mensional systems such as Hubbard-type models (in a 
loose sense) on chains and ladders. Many of these is- 
sues, but not all, can be studied in quasi-one-dimensional 
systems. However, not all of these questions can be ad- 
dressed in one dimension as the physics may be quite 
different. For instance the two-dimensional spin liq- 
uid states in two dimensions have very specific features 
with no coun terp art in one dimension (not even in lad- 
ders) [23|, EL I23 . Likewise, the description of a doped 
one-dimensional Mott insulator at weak coupling is a 
Luttinger liquid while at strong coupling is an incom- 
mensurate soliton crystal which is also a Luttinger liq- 
uid, albeit with strongly renormalized parameters. In 
contrast in two dimensions at weak coupling one may 
expect to find Fermi liquid pockets while at stronger 
couplings there is a host of possible liquid crystal like 
phases going from a solid to a stripe (or smectic) to a 
nematic whose behavior is markedly different from their 
one-dimensional counterparts (when they exist). Never- 
theless, and in spite of these caveats, studies of quasi-one- 
dimensional systems have yielded a wealth information 
on the physics of strongly correlated systems. 

The simplest quasi-one-dimensional system for the 
study of some the competing orders described above 
(and others) are ladder systems. Away from half filling 
Hubbard-type models on ladder systems can be reduced 
to the problem of two coupled Luttinger liquids. There 
is by now a rather extensive literature on the properties 
of coupled Luttinger liquids. These systems have been 
studied both analytically Effi EH E M I30L I31 151 |H 
HUE HI and numerically [HlillallaMla partly 
for their theoretical simplicity as well as a laboratory to 
test ideas intended to work possibly in two dimensions, 
and for their relevance to ladder compounds |43j . As it 
turns out, systems of two coupled Luttinger liquids can 
support almost all of the local orders proposed for two- 
dimensional systems and thus shed some light on them. 
It is thus interesting to investigate this setting the com- 
petition between different sorts of possible ordered states, 
to investigate their phase diagrams systematically and to 
compare with numerical results. 

In this paper we investigate the phase diagrams of two 
weakly coupled Luttinger liquids both at low doping lev- 
els and at half-filling, using bosonization and renormal- 
ization group (RG) methods. A number of authors have 
considered before many aspects of this problem (see in 
particular Refs. [H S M, HI IM |M |36| ) . Although 
many of the phases that will discuss here have been dis- 
cussed before, we also find a number of new and interest- 
ing phases as well as a number of new symmetry relations 
between some of these phases. 

One of the motivations of this paper was the recent 
suggestion that the Ising-like order parameter of the Z2 
symmetry of the DDW phase could be observed sepa- 
rately from the incommensuration associated with vary- 
ing the doping level 0, S H3 • If this was true it may be 



possible to have a stable phase on a ladder with spon- 
taneously broken Z2. Unfortunately, and in agreement 
with recent results by Fjarestad and Marston j3^, we 
find that while the DDW order parameter does contain 
an Ising-like piece (as it should) it always involves the 
charge degree of freedom which leads to incommensurate 
behavior. On a ladder this leads to correlation function 
which decays like a power of the distance. Although our 
results were derived at weak coupling we expect that this 
behavior should extend to strong coupling as well (with 
the usual large but finite renormalizations of velocities 
and exponents.) However, in two dimensions this implies 
at least two (and possibly more) possible and distinct 
phases: a Fermi liquid like DDW phase with pockets 
0, 0j an d a smectic (or stripe) phase with DDW or- 
der. We also find that it is quite hard to reach this phase 
in a ladder system, at least within a naive derivation of 
the effective bosonized theory from Hubbard-like micro- 
scopic models, which we summarize in the Appendix A. 
Recent, unpublished, numerical simulations by Troyer, 
Chakravarty and Schollwock have reached similar 
conclusions although in a regime where the couplings are 
larger. These authors find exponentially decaying corre- 
lations and hence only short range DDW order, which 
means that the simulations reflect a quantum disordered 
phase (of the type described below). (See also the recent 
work of Stanescu and Phillips |3.) 

The inter-twinning of charge order with some other 
sort of order (with a discrete symmetry group) is ob- 
viously not peculiar to DDW order. This is a rather 
generic situation which leads to interesting phases. It 
also happens for instance, and this is well known, to the 
Spin-Peierls or dimerized phase which upon doping in two 
dimensions it also becomes either a Fermi liquid driven 
by Fermi surface pockets at weak coupling, or a liquid 
crystal phase, such as a stripe state, at intermediate and 
strong coupling. One such example is a bond-centered 
stripe state which was considered at some length by Vo- 
jta, Zhang and Sachdev [2l|, or a site centered stripe of 
the type considered by Granath and co-workers £| which 
has a rich phase diagram. In a ladder system these phases 
are Luttinger liquid which cannot be qualitatively distin- 
guished from their weak coupling counterparts. 

We also find a number of interesting symmetries 
relating pairs of these phases. We find that, away 
from half filling, the charge density wave phase (CDW) 
with the spin-Pcircls phase (SP) (or bond-density wave 
(BW)), and a new diagonal current phase (DC) (de- 
scribed below) with the commensurate DDW phase, form 
two doublets under the continuous symmetry of slid- 
ing the charge profile, represented by the uniform chi- 
ral shift of the charge Luttinger field <fi c +: <fi c + — > 
4> c + + a (mod 4> c ± — > — 4>c± (where the real num- 

ber a is an arbitrary phase), i. e. a chiral translation 
on a circle and a reflection. This continuous symmetry 
group is non-Abelian and it may be denoted by C CX1V , 
in Schoenflies' symbols. Since in one dimensional quan- 
tum systems continuous symmetries cannot be broken 



2 



spontaneously, they can only exhibit at most quasi long 
range fluctuating order and power-law correlations. How- 
ever, at half-filling Umklapp terms break the continu- 
ous symmetry Coq V down to the finite group C^ v , i. e. 
4> c + — ► 4> c + + n^/n (mod and <p c ± — ► — <fi c ±. Hence 

at half filling these symmetries can be broken sponta- 
neously leading to true long range ordered (LRO) Ising 
type phases. In addition we also find four quantum disor- 
dered Haldane-like phases whose low energy physics can 
be described by a suitable 0(3) non-linear a model. In 
these phases there is a spin gap which remains present 
away from half filling. In this regime these phases are 
Luther-Emery liquids. There is considerable numerical 
and analytic evidence for these spin-gap phases which in 
agreement with our conclusions [33, l38l 13^ . |43| . We 
also discuss in detail the nature of the quantum phase 
transitions found at half-filling. 

This paper is organized as follows. In Section ITT1 we 
present the effective Hamiltonians and the order param- 
eters used below to characterize the different phases in 
their bosonized form. In Section IlIII we use a renormal- 
ization group analysis and the known strong coupling 
behaviors of the effective theory at low doping level to 
construct a phase diagram. In Section Hvl we do the same 
type of analysis in Section ITTT1 but at half filling. In Sec- 
tion [V] we present our conclusions. In Appendix [S] we 
relate the parameters of the effective bosonized theory 
with those of the extended Hubbard model on the lad- 
der, and in Appendix 151 we give explicit expressions for 
the order parameters of interest in terms of the bosonic 
fields. 



II. MODEL HAMILTONIANS AND ORDER 
PARAMETERS 

We begin with two coupled one-dimensional chains. To 
a large extent we will follow the approach used by Schulz 
in Ref. [2jJ. We consider first the non- interacting limit, 
and diagonalize the kinetic part in terms of "bonding" 
and "anti-bonding" bands (denoted by 1 and 2 respec- 
tively), i. e. symmetric and antisymmetric under the 
exchange of the chain labels. Including nearest neighbor 
(NN) hopping, the non-interacting dispersion relations 
are just ei a (k) = — 2tcosk =F t±(i = 1,2), where t± is 
inter-chain hopping integral. This approach makes sense 
if t± is large compared to any of the dynamically gener- 
ated gaps of the system, i. e. in the weakly interacting 
limit. 

To first order in the doping level 5, the Fermi wavevec- 
tors of two bands are, respectively, kfi^a = 7r(l — <5)/2± 
sm~ 1 (t_i/2t), and the corresponding bare Fermi velocities 
are w/i^/o = \J 4i 2 — t\ ± t±5n/2, where a is the lattice 
constant which will serve as the short distance cutoff in 
the bosonized theory. We will consider the regimes of 
both low doping and half-filling (discussed in Section ITTT1 
and Section llVI respectively) and we will assume that t± 
is not necessarily small. At half filling where the Umk- 



lapp processes dominate, the system has the particle-hole 
symmetry 



Vfi = Vf2 and kji+k 



n 



(2.1) 



Away from half-filling, we will assume that the doping 
level S is large enough to suppress the effects of all Umk- 
lapp processes (See section llTHl . However, if S is relatively 
small, the relation Eq. 12. II still holds approximately. In 
this regime the difference in their Fermi velocities does 
not play a very important role (see however the discus- 
sion in Ref. [45| )• However as the filling factor of one 
of the bands approaches zero, the respective Fermi ve- 
locity becomes very small and the physics is somewhat 
changed. In this limit there is an enhancement of the 
processes leading to the formation of a spin gap 0, • 
Since we will also find spin gap phases we will ignore 
here this special regime since it leads to the same physics 
(albeit with very different parameters). 

The effective theory consists then of two coupled Lut- 
tinger liquids, for the bonding and anti-bonding bands, 
and a set of perturbations, which we describe below, each 
associated with a particular coupling constant. In Ap- 
pendix^ we will relate these coupling constants with the 
interaction parameters of an extended Hubbard model 
on a ladder with hopping amplitudes t and on-site 
Hubbard repulsion U, and Coulomb interactions Vii (on 
the chains), V± (on the rungs) and Va (along the diago- 
nals of the elementary plaquette), as well as the exchange 
Heisenberg interactions Jy (on the chains) and J± (on the 
rungs). 

We bosonize the effective theory by introducing a 
charge bose field and a spin bose field for both the bond- 
ing and anti-bonding Fermi fields, <t> v %, where i = 1,2 
and v = c, s, where c and s label charge and spin modes 
respectively. These fields are mixed under the effects 
of various interactions, in particular the backscattering 
coupling of the respective charge and spin currents and 
densities. The bosonized theory is diagonalized in terms 
of the even and odd combinations of bose fields from each 
band: <p v ± = ± 4>v,2)/V2, U ± = {9 v ,x 
V = c, s. 

The quadratic parts of the Hamiltonian density has the 
standard "universal" form: 



T~(-c,± = 



7~L S ± — 



Vc,± 

2 

V s ,± 



Kc,±K. ± 



1 



1 



K. 



s.± 



(dx4>s,±) 2 



(2.2) 



where Tl^± are the momenta canonically conjugate to 
the bose fields 4> u ,±- The effective Luttinger parameters 
and velocities v c ,± and u S) ± are given by 



K, 



c± 



V c ,± 



1 2nv f =p 9c± 
2-kvj ± g c ± ' 



[ 2tt ' 



K, 



s± 



' 2nv f ± g s ± 
2nvf =F 9s± 



v s ,± 



- . „2 



[ 2tt ] ' 



(2.3) 
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terms appear as 



A. CDW 



B. SP 



CDC 



D. DDW 



FIG. 1: Four Ising type phases. A. charge density 
wave(CDW), B. spin-Peierls (SP); C. diagonal current (DC), 
D. d-density wave (DDW). Their triplet analogs are denoted 
as SDW, SP'; DC', DDW' respectively. 



where Vf = (vfi + Vf2)/2, The coupling constants 
9c± , 5s± correspond to forward-scattering non-chiral cou- 
plings of the charge and spin currents, and are already 
taken into account in the quadratic terms. Here we have 
ignored the effects of straightforward effects of forward- 
scattering chiral couplings, since they only renormalizc 
Fermi velocities, and modify the naively determined val- 
ues of the Luttinger parameters. Also notice that these 
expressions can only be taken seriously at weak coupling. 
At intermediate and strong couplings there are also finite 
but significant renormalization of both the Luttinger pa- 
rameters and of the velocities. 

Let us now discuss the non-quadratic, interaction 
terms. Throughout we will use Major ana Klein fac- 
tors obeying the convention 77j (1)77^ (1)7772)771 (2) = 1. 
The backscattering and pair tunneling terms yield the 
bosonized expressions 



Hi 



+ 



COS V^4 l s+ 



2(na) 2 
cos y/4n9 c 



2(ira) 



(gi cos v Anips^ — <?2 cos v 47r(9 s _ 
(53 cos V4tt9 s - + gi cos V^n<j) s ~ 



+g 5 cosV47T0 s+ ), 



(2.4) 



where 9 Ct ± and 9 Sy ± are the dual fields of the charge 
bosons 4> c ,± and spin bosons (j)s,± respectively. Terms 
labeled by the effective coupling constants g% and gi orig- 
inate from the intra-band and inter-band back-scattering 
interactions —gi{JiR Jil + 1 — » 2) and —92(Jih^2L + 
1 «-» 2) respectively. The terms labeled by the couplings 
53, (74 and t/5 represent singlet and triplet pair-tunneling 
processes X s (A\A 2 + h.c.) and A t (A{A 2 + h. c. ) with 
<73 = 2At, .94 = A s + At and 175 = X s — At. Three con- 
ditions, required by the SU (2) spin rotation invariance, 
relate the spin current and triplet tunneling couplings: 
g s ± = (51 ± 52V2 (see also Eq. and g 5 = g 4 - 53- 

Near half-filling, the following additional Umklapp 



TL U 



cos(v 47T0 C+ — 28nx) 



(g uc cos V 47t6», 



2(7ra) 2 

g u3 cos \/47r6> s _ - g ui cos V^(t>. 
g u5 cos V47r0 s+ ) 



(2.5) 



The term with coupling constant g uc is the so-called "77 
pair" tunneling processes, i. e. tunneling of Cooper pairs 
with momentum 2k /, which has the form m^mia +(1 — > 
2)+h.c, where to^l — iPr,l^R,li- The terms with cou- 
pling constants g u3 , g U 4 and g u ^ represent the couplings 
between the respective CDW and spin density wave 
(SDW) couplings on each chain: X cdw (N^ (1)N^ (2) + 
h.c), X sdw (N^(l)N^(2) + h.c), where N(i) is the 2k F 
CDW order parameter of chain 7 = 1,2, and N(i) is the 
2Uf (Neel) SDW order parameter of chain 7 = 1,2. The 
coupling constants are g u3 = ~X sdw ,g U 4:,5 = (2X cdw qp 
Xsdw/2). Due to the SU(2) spin symmetry the condition 
9u5 = gui - g U 3 also holds. 

For the two- leg ladder, we only consider where repul- 
sive interactions dominate, which implies that the bare 
values of the effective Luttinger parameters are in the 
regime K C JQ) < 1, if c _(0), # a _(0) ~ 1. Compared 
with Ref. I23, K c -(0), K s ^(0) are not necessarily 1, for 
here they are determined by off-site interactions (see Ap- 
pendix EJ . 

Bosonic expressions for various order parameters are 
given in Appendix [5] I n the particle-hole (p-h) chan- 
nel, the possible singlet fermionic bilinear forms which 
break the translational symmetry are the order param- 
eters for the CDW and SP, DC and DDW operators as 
shown in Fig. ^ The CDW and SP order parameters are 
proportional to the real and imaginary parts of the sym- 
metric bilinear ip\ La ip2Ro + i^La^^R" ^ whereas the DC 
and DDW order parameters are the real and imaginary 
parts of the anti-symmetric version of this bilinear. 

From their bosonic representations, we find that all 
four order parameters transform non-trivially under the 
symmetries broken in their associated phases (or ground 
states). Thus, for instance the SP and DDW order 
parameters are odd under the Z2 symmetries broken 
spontaneously by the SP and DDW phases. However, 
in all four cases, these order parameters also involve a 
phase factor (or vertex operator) of the charge boson 
4> c ,+- Hence these order parameters also transform non- 
trivially under shifts of the charge boson 4> c .+, i- e. uni- 
form displacements of the charge profile. This depen- 
dence means that the discrete symmetries, broken spon- 
taneously in these phases with long range order, are inter- 
twinned with the continuous symmetry of the incommen- 
surate doped state. Consequently, these order parame- 
ters do not truly acquire an expectation value but in- 
stead only display power law correlations. Also, while 
it is possible to write down bosonic expressions for op- 
erators which transform only under the discrete symme- 
tries broken by these phases, their fermionic versions are 
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strongly non-local. Hence, we conclude that these orders 
are always incommensurate. 

We also find that these order parameters also form two 
doublets of the Coo V group. Similarly, their triplet coun- 
terparts SDW, SP*, DC*, DDW* are proportional to real 



and imaginary parts of ip[ La (% ) a p'>l>2Rp±'iPl La (% )api>inp 
respectively (where the label t means triplet). In the 
particle-particle (p-p) channel, the s and d-wave pair- 
ing order parameters are A s ^ d = '(^lLa-^im ± 
ip2L<r' l PiRa) • I n the next section we identify the stable 
fixed points of the renormalization group (RG) flows for 
the these phases associated with these order parameters. 

Some of the order parameters discussed above have 
been investigated before in Ref. (25^ |3(j although under 
different names. For example, our CDW, DDW, SDW, 
SSC and DSC order parameters are called CDW"*, 
OAF, SDW*, SC S and SC d there. We note that in 
a recent paper Ref. ^(|, the phases that we label as 
DDW, SP and DC are called called SF, P-density wave 
(PDW), F-density wave(FDW) respectively. 

Finally, in Eq. 12.11 we ignored the effects of the fol- 
lowing terms 

ah c = (Av f + ^)d x <j> c+ d x <j> c - + (Av f - ^)n c+ n c _, 

7T 7T 

(2.6) 

ah s = (Av f - — ;/v.WVx + (Av f + ^)n s+ n s _, 

7T 7T 



2(7ra) 2 



sin v A~K(f) s - sin v 47T(/>s+ 



(2.7) 



ATLum — 



sin(V47r0 c+ — 2Sttx) 



(Ag uc cos v 4:Tt6 c - 



2(7ra) 2 

— Ag U 3 cos v47r0 s _ — A<7„4 cos \fkK<$> s - 

- A.g u5 cos\/47r0 s+ ), (2.8) 

where Avf = 5nt±/2 and all other residue coupling con- 
stants varnish linearly with doping near half-filling as 
given in Appendix A. The quadratic residual terms in 
Eq. 12.61 and Eq. 12.71 are marginal perturbations, and 
they slightly change the scaling dimensions of various 
operators in Eq. 12.41 12.51 Because they are small, we 
do not expect that they can change the stable RG fixed 
points associated with various phases qualitatively. For 
the term of Ag s in Eq. 12.71 6 S + is fixed around or 
y / 7r/2 at all the stable fixed points (see Table HI and IIIII 
below). The residual Umklapp terms in Eq. 12. 81 are irrel- 
evant away from half-filling. At half-filling, 9 C+ is fixed 
at 0F/2 (see Table UnJ). Thus we conclude that all the 
non-quadratic operators are irrelevant at all the stable 
fixed points. Balents and Fisher used a perturba- 
tive RG of the fermionic theory and found that a spin 
gap phase develops near half-filling, which is consistent 
with the argument given above. On the other hand, the 
continuous C voo symmetry is preserved away from half- 
filling where the Umklapp terms are irrelevant. Thus the 
conclusion that the CDW and SP, DDW and DC order 



parameters are incommensurate and thus exhibit quasi 
long range order is not affected by these terms. How- 
ever, these residual terms do affect the boundaries among 
phases. 



III. PHASE DIAGRAM IN THE 
INCOMMENSURATE REGIME 

We will now investigate the phase diagram in the in- 
commensurate regime, but only at low doping. In this 
regime, the Umklapp processes are cut off at a high en- 
ergy scale of 2irvf5/a, and can only yield renormaliza- 
tion of the parameters, such as the velocities, coupling 
constants and Luttinger parameters of the low energy 
effective theory. The contributions from the Umklapp 
terms in the RG equations away from half-filling ^7j are 
given in terms of Bcsscl functions, which oscillate when 
a energy scale lower than that of the Umklapp process 
is reached. At this scale, the effects of these terms can 
be neglected. Below, we begin directly at the low en- 
ergy scale with all the coupling constants and Luttinger 
parameters already renormalized by the Umklapp terms. 

We will investigate the role of the remaining interac- 
tions by means of a one-loop renormalization group (RG) 
analysis combined with semiclassical arguments. In this 
regime, the charge boson <p c .+ essentially decouples and 
remains gapless. Thus, to one loop order, the Luttinger 
parameter A Cj+ does not flow. (This argument is not 
completely correct: there are always irrelevant couplings 
which do lead to finite renormalizations of AT C:+ ; these 
effects do not show up at one-loop order.) 

The one-loop RG equations for the coupling constants 
gi through g$ and Luttinger parameters K c _ and K Si ± 
are 



dl 




-^(gi + gt + gl) 


dK s+ 
dl 




K 2 

8jT 2 \9i + .92 + 


dK s _ 
dl 




K 2 1 

v s- / 2 1 2\ 1 1 / 2 

g7r2 C9i+.9 4 )+ 87r2 (.92 


dgi 
dl 


(2 


K s+ 3 f 


dg 2 
dl 


(2 


• A s > + 2n 


dg3 
dl 


(2 


1 1 x , 5235 

K c - K s - )9 * 2n 


dg 4 
dl 


(2 


1 v \ digs 

A c _ Z7T 


dg5 
dl 


(2 


AV Ks+)95 2, 



2tt ' 

(3.1) 

where / = ln(L/a) with the length scale L. 

Along the S , J7(2)-invariant manifold for the spin cur- 
rent and pair tunneling terms, the RG equations can be 
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simplified to 






dK c _ 
dl 




-gl + {93-6 


dg s + 
dl 




2 \ (93 

+ 91-) 


dg s - 
dl 


1 

--gs+g s 


2 2 
477 47T 




1 




AV + 




1H- 


1 




AV + 


2n 



-in 
)93 - 

)94 - 



(g s + - g s -)g4 

(g s + + 9s-)93 



2vr 



(3.2) 



with 



9s+~ 



37(93-54+35) = (l-i^ + ^)(93-94+<?5) = (3.3) 



dl 



K, 



2n 



These equations are invariant under transformations 

(5i j 92, 93, 94) -> (91,92,-93,-94) (92,9i,94,93)- 
This means phase boundaries must also have such sym- 
metries. 





9i,32 


33,54,35 


4>s+ 


4>s- 


6> s - 


Order 


dimension 


1 


0, — 00 


+00, 0, — 00 





1 


2 


CDW+SP 


A c+ /4 


2 


0, —00 


—00, 0, +00 


2 


1 





DC+DDW 




3 


—00, 


0, +00, +00 


2 


2 


/ 


DSC 


l/(4^c+) 


4 


—00, 


0, — 00 , —00 








/ 


SSC 


V(4^ c+ ) 



TABLE I: Stable fixed points and corresponding quasi long 
range orders away from half- filling, with {0 C -} = and g$ = 
9i — 93 (required by SU(2) invariance.) 





9s+,gs- 


33,34,35 


<j>s + 


Transition 


1 


0,0 


+00, +00, 


unfixed 


CDW+SP <-> DSC 


2 


0,0 


—00, —co, 


unfixed 


DDW+DC <-> SSC 


3 


— co, 


— CO, +co, +00 


0F/2 


DDW+DC <-> DSC 


4 


—co, 


+00, —co, -co 





CDW+SP <-> DSC 



TABLE II: Critical phase boundaries and unstable fixed 
points away from half-filling, also with (9 C -) = and g$ — 
34 - 33- 



For "bare values" of the Luttinger parameter K c _ (0) ~ 
1, the marginally relevant RG flow of Eq. (|3.2|l is such 
that a gap develops in the c— sector, which scales like 
to c _ sa exp(— l/g(0)), where g is the most relevant one 
among the marginally relevant perturbations 33, 34 and 
35. In this regime K c — flows to large values and, thus 
from now on we will set 1/K C - = 0. In this phase the 
operator cos(V4~7r# c ,_) acquires a non- vanishing expec- 
tation value which classically are just ±1. Hence in this 
phase the dual field takes the values 9 C - = 0, \/n/2, 
which are related to each other by a Z2 symmetry |3fij . 
In what follows in this section, we will choose the value 
(0c-) = 0. 



From now on we will use the set (31,32,93,94) to rep- 
resent the stable fixed points of Eq. (|3.2|) , which are sum- 
marized in Table [i] At the fixed points (0, —00, +00, 0), 
the inter-band back-scattering coupling constant 32 is 
relevant, while the intra-band back-scattering coupling 
constant 31 is irrelevant. Both X s and At are relevant 
and satisfy the relation X s — —At. By direct inspec- 
tion of their scaling dimensions, we find that X s and At 
are more relevant than g^- The resulting phase depends 
on where the RG flows go. When g 3 — > —00, the ex- 
pectation value of (f> s + and 6 S _ asymptotically take the 
values ((j> s +) = V^/2 and (9 S -) = 0. This is the sta- 
ble fixed point for either the DDW phase or the DC 
phase. However, this is true only quasi-long range or- 
der (QLRO) due to the strong fluctuations of the gapless 
charge boson <p c .+- in this phase these order parame- 
ters have scaling dimension K c ^ + /A. Conversely, when 

33 — > +00, (<j) s +) = 0, and (9 S -) = \/n/2. Hence, at this 
fixed point we would have (naively) either a CDW phase 
or a spin-Peierls (SP (or dimerized) phase. Here too there 
is only QLRO and the associated order parameters also 
have scaling dimension K c + /A. 

We conclude, in agreement with the recent results of 
Ref. j3|| , that the chiral translation symmetry in the field 
4> c ,+ i m other terms due to the charge incommensurabil- 
ity, there is no true long range order of the DDW or- 
der but only (incommensurate) power law correlations. 
We further can see that the DDW and DC phases ( and 
also the CDW and SP phases) form doublet represen- 
tation under the Coo V group and are thus degenerate. 
Equivalently, the DDW and DC order parameters can 
be regarded as the real and imaginary parts of a single 
complex order parameter which can thus be rotated con- 
tinuously into each other. The same relationship holds 
for the CDW and spin-Peierls order parameters. Thus 
both stable phases, CDW+SP and DDW+DC, have a 
continuous U(l) symmetry. Naturally, since the ladder 
is a one-dimensional system, this symmetry is no truly 
spontaneously broken as there are only power-law cor- 
relations for these order parameters. However, we will 
see in Section llVl that at half filling, the Umklapp terms 
break this symmetry explicitly from U(X) down to Z2 
leading to additional Ising-like phase transitions. 

Similarly, we also find that A s is more relevant than g\ 
at (—00, 0,0, ±00) while 32 and At are irrelevant. When 

34 — > +00, {4>a+) ancl (0a-) are fixed at \pnj2. Thus, 
d-wave superconductivity (DSC) is the leading QLRO 
and its order parameter has scaling dimension 1/(AK C+ ). 
Conversely, when 54 — ■> —00, (4> s +) and (<f>s-) are 
fixed at 0, s-wave superconductivity (SSC) is the leading 
QLRO and its order parameter also has scaling dimen- 
sion 1/(4A' C+ )- 

Let us consider now the phase boundaries and the 
nature of the phase transitions between these possible 
states, at 3 S _(0) = 0. In this regime it is more natu- 
ral to represent instead the unstable fixed points with 
{g s +, 3s- , 33, 34)- The RG flows starting with 3 s +(0) > 
0,3 S _(0) = and 33(0) = 34(0) = g > evolve towards 
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0.4 







VX\ DSC 




DDW+DC Xi 


CDW+SP 


ssc 









-0.4 -0.2 0.2 0.4 



93 

FIG. 2: Phase boundaries with positive initial value of g s + 
( g s +(0) — 0.2) and different initial values of g s -(0) with 
dashed line (g s _(0) = 0), solid points (g s _(0) = 0.1) and 
triangles (g 3 _ (0) = -0.1). Phase boundaries of CDW+SP <-> 
DSC, DDW+DC <-> SSC become of first order for g s+ (0) < 0. 



the fixed point at (0, 0, +00, +00). Here, the field <f) s+ 
becomes free, K s ± — > 1, and the residual interactions 
reduce to 

nl es = n , 9 , 9 (cosv^6> c _)(cos\/47rfl s - +cos%/4^0 s _), 

(3.4) 

where g* means the renormalized value of g. At this 
fixed point K s - — ► 1, and both perturbations are di- 
mension 1 operators. This system is invariant under the 
duality transformation S _ <-> # s _. This model has been 
studied extensively in the literature |29t l48| . It is equiv- 
alent to a theory of two Ising models. If the coupling 
constant in front of both operators is the same, as it 
is the case in Eq. I|3.4[) . one of the Ising models is at 
its critical point. Equivalcntly, it can be regarded as a 
theory of two Major ana fcrmions, one of which is mas- 
sive. Hence this fixed point is in the universality class 
of the two-dimensional classical Ising model. The Ising 
order and disorder operators are given by sm.y/Tr^) s _ and 
sm-^fnOs- respectively. At this fixed point, both oper- 
ators have scaling dimension 1/8, as they should at an 
Ising transition. A small perturbation making (73 > 174 
or 93 ^5 94 causes a flow towards the CDW+SP or DSC 
fixed points respectively. Thus 53 = 54 > is the phase 
boundary between the phase CDW+SP and a d-wave su- 
perconductor at g s - (0) = and g s +(0) > 0. 

However, if the RG flows begin with g s +(0) < along 
this direction, then the field S _ is no longer critical. Ac- 
cording to Eq. 13.21 in this regime g s + is marginally rele- 
vant and g s+ — > —00, with (73 = g± > and g\ = g2 < 0. 
At this fixed point, the fields 9 C _ and <j) s+ acquire non- 
vanishing expectation values, and the residual interac- 




FIG. 3: RG flows in the 3D parameter space with g s +(0) > 0. 
The dashed lines mark the critical surface. \g^\ wins over \g&\ 
on the left of the surface, and |<j4| wins over \gs\ on the right. 
On the critical surface, the RG trajectories flow to the line 

\gs\ = M- 

tions at this fixed point reduce to 

= COb s — [g* (cos V4^6 C -) + g{ (cos \fl^4> s +)) 

+ C ° 2(^)2^ (ff3(cos\Z4^6> e _) - g ;(cos V4^> s+ )). 

(3.5) 

At this stage of RG, the renormalized couplings sat- 
isfy g\ = g\ and g\ — g\. Once again we can take 
(cos(\/47r# c _)) = 1 (the renormalization of its amplitude 
can be absorbed in a redefined coupling constant). This 
effective theory has the same form as Eq. I|3.4|) . Hence 
this is also a theory of two Ising models. However, un- 
like Eq. H3.4fl the amplitudes of the two dimension one 
operators are not equal. Hence, generically both Ising 
models are off-criticallity (or equivalently both species 
of Majorana fermions are massive). This corresponds to 
a finite correlation length and finite energy gap at the 
phase boundary. Hence, in general this is a first order 
transition. If (4> s +) = 0, then the term of # s _ wins over 
that of 4> s - and g s - — ► +00 in the next step RG trans- 
formation. Conversely if (4> s +) = then the term 
of 4> s - wins over that of # s _ and g s - — > —00 in the next 
step RG transformation. Finally, RG flows evolve to the 
CDW+SP fixed point in the former case while towards 
the DSC fixed point in the latter. Thus for 33 = 54 > 
and g s+ (0) < 0, the phase transition at the boundary of 
CDW+SP <-> DSC becomes first order as the correlation 
length is now finite. However, a second order transition 
is also possible here too. If the spin boson cf> s+ is quan- 
tum disordered, then (cos(v / 47r</ ) s+)) = and once again 
we get an Ising critical point of the same kind discussed 
above. Hence the general conclusion is that this phase 



7 



boundary may be at a second order transition (with Ising 
criticallity) or at a first order transition, with an Ising- 
like tricritical point in between. Similarly, g 3 = 54 < 
at ffs „ = is the boundary of DDW+DC <-> SSC, which 
is critical and leads to the fixed point at (0, 0, — oo, — oo) 
or to a first order when g s+ > or g s+ < respectively. 

Another pair of fixed points (— oo, 0, +oo, ±oo) control 
the phase boundaries of the DDW+DC «-> DSC tran- 
sition at 33 = —54 < 0, where (4> s +) — 0, and at the 
phase boundary of the CDW+SP <-> SSC transition at 
.93 = -.94 > 0, where (cf> a+ ) = y/n/2; g s+ -> -oo no 
matter what its initial value is. The residual interaction 
for the s— sector is still described by Eq. (|3.5II but now 
with 53 = —g%. Thus, the amplitudes of cosV^tt^s- 
and cos >/4n9 s - are kept equal and this phase bound- 
ary is also in the universality class of the Ising critical 
point. The Ising order and disorder operators can be de- 
termined accordingly. The critical phase boundaries are 
summarized in Table ITT1 

The initial value g s - (0) has important effects on phase 
boundaries. In Fig.|2we present the result of a numerical 
integration of Eq. I|3.2[) for <? s +(0) > 0; <? s _ > favors the 
growth of |<73| but disfavors that of |<?4|, and conversely 
g s - < favors the growth of I34 1 but disfavors that of 
|.g3 1 . Let us begin with the case g s - (0) > 0. For 1 1?3 (0) | < 
|54(0)|, at first <7 S _ decreases, then it reaches a positive 
minimum and finally it increases. Thus \gz\ increases 
faster than \g 4 \ and eventually it wins over it. However, 
if 1.93(0)1 -C 154(0)1, g s - decreases monotonically to nega- 
tive values and I54I still wins over \g 3 \. As a result, both 
regions of the phase diagram with DDW+DC order and 
CDW+SP order expand beyond the line 33 = ±34, and 
the correspondingly areas of d-wave and s-wave supercon- 
ductivity correspondingly shrink. Due to the symmetry 
of Eq. (|3.2|) . the situation is reversed for 3 S _(0) < 0. For 
an initial point located on one of these phase boundaries, 
its RG trajectory flows to the corresponding unstable 
fixed point, as shown in Fig. [3J For 3 s +(0) < 0, the effect 
of 3 S - (0) is similar, but the phase boundaries CDW+SP 
<-> DSC and DDW+DC <-> SSC are now first order tran- 
sitions and there are no accessible critical points. 

We conclude this section with some comments on the 
DDW phase which has attracted considerable interest re- 
cently. Until now there is no solid numerical evidence 
away from half- filling [Zol 14^ ] . For the two- leg ladder, we 
find (see Appendix that the DDW phase may exist 
but it is necessarily incommensurate. We also find V± 
large and positive reduces 34 and enhances g s -, which 
is favorable for the DDW phase to exist. However, a 
negative 33 with magnitude comparable to I34I is also 
needed. Thus we suggest to look for it in the regimes 
V± ~> Vd >■ Vji > 0, which has only repulsive inter- 
actions, or for VI > > Vji , which has some attrac- 
tive interactions (and thus is less physically relevant). 
These arguments agrees with the results of a recent two- 
dimensional mean-field calculation that the Hubbard 
U alone can not stabilize the DDW phase and that neg- 
ative nearest-neighbor interactions are needed. However, 



V\i < together favor d-wave superconductivity over 
the DDW state. 



IV. THE PHASE DIAGRAM AT HALF FILLING 

Let us now discuss the phase diagram at half-filling. 
The main change is the presence of Umklapp terms. 
Compared to the incommensurate case discussed in Sec- 
tion ^Qj the main difference is that at half-filling the Z2 
symmetries behind two-fold degeneracies found in away 
from half-filling are now can be spontaneously broken, 
with possible phase transitions between the CDW and 
the spin-Peierls phases, and between the DDW and the 
DC phases. Since much of the analysis is rather similar, 
here we will only sketch the main differences. 

The set of RG equations are now more complicated: 



dK c+ 

dl 
dKc- 

dl 
dK s+ 

dl 
dKs- 

dl 
dgi 
dl 

dm 

dl 
dgi, 
dl 
dgi 
dl 
dgh 
dl 

dg U c 



K 2 

c+ / 2 1 2 1 „2 



87T 



1 



8tt 2 
K1_ 



(gi + gi + gt + 9I5) 



8^5- (5? + 9l + 9I4) + (ffi 



g\ 



9I3) 



(2 
(2 
(2 
(2 
(2 



K 
K 



K s -)g x 
1 , 



5455 5«45«5 



s+ 

1 



1 



dl 
dg U 3 

dl 
dg U 4 

dl 
dg U 5 

dl 



K c - 
1 

1 

(2 - K c+ 
(2 - K c+ 
(2 - K c+ 



.92 - 

)93 



K. 
K s -)g4 

K s +)g5 
1 



2tt 


2?r 


5355 


5«35«5 


2tt 


2vr 


3255 


9us9uc 


2tt 


2tt 




g-uiguc 



3i34 , 5253 , 5«55« 



+ 



I 



)gn 



)gu3 



2tt 2tt 2tt 

. 535«3 . 545«4 . 555»5 



K s 

- K s - )g uA - 
(2 - K c+ - K s+ )g u5 - 



2vr 

929u5 

2tt 

5lg«5 

2tt 

gigm 

2tt 



2tt 

g3g uc 

2ir 

545W 
2tt 

325m3 

2tt 



2tt 



555mc 

2tt ' 
(4.1) 



We will not be interested here in solving these RG equa- 
tions in their full glory, but only in the regime where 
K c+ <C 1 and K c - ~ 1. For this range of parameters 
there are a number of useful hierarchies of scales which 
considerably simplify the analysis. 

Contrary to what happens away from half filling, the 
field 4> c+ no longer decouples due to the effects of the 
Umklapp terms of Eq. 12.51 Clearly, (f> c+ plays a role 
quite similar to that of 9 C -. Indeed, in this regime, g uc 
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is the most relevant coupling and it is associated to an 
operator with scaling dimension K c+ + 1 /K c - ■ This op- 
erator takes the RG flow close to a fixed point at which 
the field <t> c + acquires a gap approximately of the form 
m c+ w a _1 |3 nc (0)| 1/(1 " ifc + (0)) . In this regime the field 
9 C - behaves roughly in the same way as in Eq. (|3.2|l . Here 
too the coupling constant g c - flows to strong coupling, 
1/K C - — > and a gap m c _ develops in this sector as it 
does away from half filling. We will set (4> c +) = 
correspondingly (6 C -) = or y/n when g uc (0) > or < 
respectively, so that (cos \/47r<^ c _|_) ~ —(a ni c+ ) Kc +^ and 
(cos V^Oc-) « sgn( ffuc )(a m c _)V*.-(o). 

Once the fields C + and 9 C - become pinned close to 
their classical values, the effective residual interactions 
among the remaining fluctuating degrees of freedom have 
an effective Hamiltonian of the form 

COS , 1-7— , n~ n \ 

"e// = — — ^2 — (91 cos V47r</> s _ - g 2 cos V47r0 s _) 
Zyira) 

+ —. 1 , 9 (ffa cos V-hrO s - + 54 cos V-hr^s- 
2(7raj' 2 

+ g* 5 cosV^<j) s +), (4.2) 

where 53,4(0) = 53,4(0) (cos V^0 C -) - 

5«3,«4(O)(cosv / 47r0 c+ ) and g|(0) = .94(0) - 33 (0). 
If 5«c(0) is not small compared with the initial 
(or bare) values of the other coupling constants, 
this first step of the renormalization group flow is 
rather quick. In this step the marginal coupling 
constants cannot not change very much and thus 
I (cos v / 47T0c+) I S> I (cos V^r^c-)! is a good approxima- 
tion. Hence the renormalized residual couplings are 
approximately {gl,g\,gt) ~ (0*3(0), 5 u4 (0), 0*5(0)). 





Que 


91, 32 


53,34,35 


0c- 




(0.-A-) 


phase 


1 


+00 


0, — 00 


+00, 0, —00 








(/, ^) 


SP 


2 


+00 


0, — 00 


—00, 0, +00 





^F 
2 


(/ ,0) 


DDW 


3 


+00 


— 00 , 


0, +00, +00 





2 


(#>/) 


DSC+SDW 


4 


+00 


—00, 


0, —00, —00 








(0, /) 


SSC+DC* 


5 


—00 


0, —00 


—00, 0, +00 


2 


2 


(/ ,0) 


CDW 


6 


—00 


0, —00 


+00, 0, —00 


VFF 

2 





(/, ^) 


DC 


7 


—00 


— 00, 


0, —00, —00 


2 





(0,/) 


DSC+SP* 


8 


—00 


—00, 


0, +00, +00 


2 


V^F 

2 


(#,/) 


SSC+DDW' 



TABLE III: Fixed points at half-filling: Stable fixed points 
and corresponding gapped phases. We have set (4> c +) ~ 
■/tt/2. The SU(2) condition requires gi, =34—33- Phases 
1,2,5,6 have true Ising type long range order, while 3,4,7,8 are 
quantum disordered Haldane-like phases . 

The new RG equations, which control the subsequent 
RG flow, are the same as in Eq. I|3.2[) after setting 
1/K C - — ► 0. Here too the SU(2) condition g\ = g% — g^ 
is obeyed, albeit among renormalized couplings. The re- 
sulting stable phases and the phase boundaries between 
them are given in the phase diagrams of Fig. 0] and 





9uc 


9s+ 


33,34,35 


e c - 


0S + 


Transition 


1 


+00 





+00, +00, 





/ 


DSC+SDW^ SP 


2 


+00 





—00, —00, 





1 


SSC+DC* <-> DDW 


3 


+00 


—00 


—00, +00, +00 





nAF 
2 


DSC+SDW <-> DDW 


4 


+00 


—00 


+00, —00, — 00 








SSC+DC' <-> SP 


5 


— 00 





—00, —00, 


2 


/ 


DSC+SP* <-> CDW 


6 


— 00 





+00, +00, 


sAF 
2 


/ 


SSC+DDW* <-> DC 


7 


— 00 


—00 


+00, —00, —00 


& 
2 





DSC+SP* <-> DC 


8 


—00 


—00 


—00, +00, +00 


2 




SSC+DDW* <-> CDW 



TABLE IV: Fixed points at half-filling: Unstable fixed points 
which have the common fixed value g s - = 0. Here too 
(4>c+) = \A~/2, and the SU(2) condition requires g% — 34—33. 
The column on the right indicates which transition is con- 
trolled by each unstable fixed point. 

[3] The corresponding stable fixed points and values of 
pinned fields are summarized in Table ITTll The critical 
(or unstable) fixed points are given in Table ITVl Umk- 
lapp terms break the symmetry group to C4-, and thus re- 
move the degeneracy between CDW and SP phases, and 
between the DDW and DC phases. Hence, all four states 
become distinct phases with true long range order, which 
break the residual Z2 symmetry spontaneously. At the 
quantum phase transitions between CDW and SP, and 
between DDW and DC, the symmetry is U (l). 

Perturbative RG studies of Refs. |35j and 36] have de- 
scribed the CDW and DDW fixed points with the prop- 
erty that the coupling constants (written in our notation) 
satisfy 

- 02 = ±53 = ±55 = -5u3 = 0*5 = T9uc -> +00 

5i = 54 = 5u4 = 0. (4.3) 

where the upper (lower) sign holds for the CDW (DDW) 
phase. It turns out that a model with this particu- 
lar choice of coupling constants was been proposed by 
Scalapino, Zhang and Hanke (SZH) [34[ as a ladder 
model of the 50(5) theory. However, Lin, Balents and 
Fisher found that, at least to one loop order in a per- 
turbative RG |3f|, the symmetry is enlarged actually to 
50(8). Moreover these authors found, also within a per- 
turbative RG, that the 50(8) manifold is at least locally 
stable, i. e. small deviations from this trajectory con- 
verge to this trajectory under the RG flow. Interestingly, 
the 50(8) manifold is an integrable fermionic system for 
which a number of exact properties have been calculated 
using the Bethe Ansatz 49]. 50(8) is clearly a dynam- 
ical symmetry which is possible because the operators 
that are involved (back in the fermionic representation) 
are all dimension 2, they are superficially marginal but 
become marginally relevant due to fluctuations leading 
to the development of a gap. 

However, for more generic values of the coupling con- 
stants this dynamical symmetry does not necessarily 
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arise. It is not know how large the basin of attrac- 
tion of the SO (8) manifold actually is. In fact using 
bosonization methods we find that far away from the 
SO (8) manifold the scaling dimensions of these opera- 
tors begin to differ significantly from each other and thus 
evolve differently under the RG (see Eq. I|4.3jl ). In partic- 
ular, by checking their scaling dimensions, we find that 
the renormalized couplings can renormalize so differently 
from each other as 

I02I < \9a\ = \gs\ < \9u3\ = Iswl < |ffc| -> 00 
9x,9i,9ui (4.4) 

in all the four phases of CDW, DDW, SP and DC ( Re- 
call that the signs of the coupling constants change in 
some of the phases. ). Nevertheless, what is clear is that 
the spectrum found in these more anisotropic (and more 
generic) regimes is smoothly connected to the multiplets 
found in the 5*0(8) limit. In other words there is no phase 
transition separating these regimes, but the spectrum is 
organized differently. 

Let us now discuss the phase transitions between the 
CDW and SP phases, and between the DDW and DC 
and phases, and to the associated critical fixed points. 
As we noted before, these phase transitions are driven 
by the Umklapp terms, the most relevant of which is 
controlled by the coupling constant g uc . Hence, at the 
critical point separating the SP and CDW phases, and 
between the DDW abd DC phases, the Umklapp terms 
are tuned to zero. The critical fixed points coincide with 
the stable fixed points of the incommensurate CDW+SP 
phase and DDW+DC phase respectively. In both cases 
the transition is controlled by the sign of g uc . We also 
note that the renormalized coupling constant g^ has dif- 
ferent signs on both sides of this phase transition. This 
is because, close to the transition g% « 53 (cos ^/4tt9 c -), 
and (# c -) = in the SP phase while (0 C -) = m the 

CDW phase. The same is true for the phase transition 
between the DDW and DC phases. 

It can be shown that, if only charge interactions are 
considered |50j . then g uc — g u ^ at the bare level. In this 
regime the CDW and SP phases are more easily accessible 
than the DDW and DC phases. There is strong numerical 
evidence for a commensurate DDW phase at half-filling 
in a t-J-Hubbard ladder |4lJ who included Heisenberg- 
like exchange interactions at the microscopic level. It is 
easy to see that although the inclusion of microscopic ex- 
change interactions does not lead to a different low energy 
theory, it changes the strengths of the different effective 
couplings. In particular it makes the DDW phase more 
accessible. For simplicity, we discuss the conditions of 
the commensurate DDW phase on the SZH ladder which 
only includes nonzero interactions U, V± , Jj_ . The cou- 
pling constants are given in the weak interaction limit 
in Appendix A. Let us suppose that V± and J± > 0. 
First of all, we need positive g c+ to set up the overall 
repulsive interaction i.e., U + 2V± > 0. A large J± helps 
to make g uc > and g u z < simultaneously i. e. , 
j,J± > U — V± > — I Jj_. But J± can not be too large, 



otherwise negative g s _ suppresses the DDW phase. For 
\g U 3\ > which can be achieved with U < 0, this 

phase is stabilized. But |?7| can not be too large, oth- 
erwise g c + would become negative. The re gion where 
the commensurate DDW was found in Ref. 4lJ] agrees 
with this analysis. Again, we need to keep in mind that 
this naive analysis only makes sense in the weak coupling 
limit, which also neglects effects from many irrelevant op- 
erators. Thus, we do not expect this analysis to give a 
precise location of the phase boundary. 

Now we discuss the remaining phases and phase tran- 
sitions. Upon a careful study of which fields become 
pinned and what are their allowed expectation values, 
we conclude that the remaining four phases are actually 
quantum disordered Haldane-like states. For example, 
there is a phase in which d-wave superconductivity and 
the SDW order parameters (DSC+SDW) are quantum 
disordered. The order parameter for DSC is very sensi- 
tive to fluctuations in the c+ sector since Odsq oc c + 
Similarly, the x, y and z components of the SDW order 
parameter are controlled by fluctuations in the s± sector 
since O sdw oc (sin( v / 7r6' s _), sin(^ s+ ), cos(y/Tr9 s+ )). 
At this fixed point the fields 6 C+ and 9 S ± are not pinned 
and fluctuate wildly. Nevertheless, the remaining fields 
in the expressions for these order parameters do pro- 
vide for a finite amplitude even though the fluctuations 
of both phase and orientation are so strong that the 
system is quantum disordered. The analysis of other 
three phases, s-wave superconductivity and triplet DC 
(SSC+DC*), d-wave superconductor and triplet spin- 
Peierls (DSC+SP*), and s-wave superconductor and 
triplet d-density wave (SSC+DDW*), is similar. Because 
of large charge gaps, the low energy physics of their spin 
sector, may be described by the corresponding 0(3) non- 
linear a model without a Berry phase term, which is 
quantum disordered. 

The phase transition between the DSC+SDW phase 
and the DSC+SP* phase (see Fig. H and Fig. 0) is 
the commensurate limit of the d-wave superconductor 
found away from half filling. A similar relation holds 
for the phase transition between the SSC+DC* phase, 
the SSC+DDW* phase and the s-wave superconductor. 

Finally, let us discuss the unstable fixed points with 
I #3 1 = Isll ~~ * °°i5s- = 0, summarized in Table HVI 
The RG flows starting from the phase boundaries with 
g s +(0) > evolve towards these fixed points. At these 
phase boundaries, the order parameters for CDW, SP, 
DC and DDW have power-law correlations and have scal- 
ing dimension 3/8 at the fixed points denoted by 1,2,5, 
and 6, and scaling dimension 1/8 at the fixed points de- 
noted by 3, 4, 7 and 8 (see Figs. 0] and |5J). On these 
phase boundaries the d-wave and s-wave superconduct- 
ing order parameters are quantum disordered. Simi- 
larly, the SDW, SP*, DC* and DDW* order parameters 
have power-law correlations and their scaling dimension 
is 3/8 at the points 1,2,5,6 but are quantum disordered 
at points 3,4,7,8. For g s -(0) — 0, at these phase bound- 
aries the renormalized couplings satisfy | g% \ = \g%\ as 
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FIG. 4: Stable phases and phase boundaries at half-filling 
with g s _(0) = 0,g s +(0) > 0, and g U c(0) > 0. Phase 
boundaries 1,2,5,6 are represent first order transitions when 
g s +(0) < 0. The critical fixed points for the transitions from 
phases in this figure to their counterparts in Fig. |S] are anal- 
ogous to those of Fig. [5] 



B 7 §4 5 




6 8 



FIG. 5: Phase diagram and boundaries at half-filling with 
gs_(0) = 0, g s +(0) > and g U c(Q) < 0. Notice that here we 
use — <?•; and —g% as the x,y axes. 

before. Nonzero g s _ (0) also has similar effects on these 
phase boundaries: g s -(0) > (<)0 favors phases CDW, 
SP, DC and DDW (DSC+SDW, SSC+DC*, DSC+SP* 
and SSC+DDW* ) respectively. When 5s+ (0) < 0, the 
situation is similar except that transitions 1,2,5,6 become 
the first order and there are no corresponding unstable 
fixed points. 

V. CONCLUSIONS 

In summary, in this paper we studied the problem of 
competing orders in two- leg ladders, which were mapped 
to two-coupled Luttinger liquids with p-h symmetry at 



both low doping and at half-filling. We used (Abelian) 
bosonization and RG methods to study the phase dia- 
grams of these ladders both at half filling and at low 
doping. Stable and unstable fixed points of the RG 
flows with the corresponding phases and phase bound- 
aries were investigated in detail. First order transitions 
when g s +{0) < are found and the effects of g s -{0) 
to phase boundaries are discussed. The Coo V symmetry 
makes CDW and spin-Peierls, DC and DDW degenerate. 
In the absence of Umklapp terms there is incommensu- 
rate quasi long range order. These degeneracies are re- 
moved at half- filling where true long range order appears. 
Power law fluctuating c?-wave and s-wave supreconduct- 
ing phases at low doping levels become quantum disor- 
dered at half-filling, with finite amplitudes among DSC, 
SSC and SDW, DC*, SP*, DDW* respectively. Sugges- 
tions on how to best find these phases in numerical sim- 
ulations were given. 

After this paper was submitted for publication, we be- 
came aware of the work by Tsuchiizu and Furusaki on a 
very similar model (at half-filling) 46] . In this work these 
authors also obtained the same eight insulating phases we 
found here at half-filling. Also after this work was sub- 
mitted, we learned of the numerical work by Schollwock 
et. at. 51] on a DMRG study of a similar ladder model 
away from the half filling. At low doping these authors 
found that their results are consistent with an inhomo- 
geneous picture of the doped state in which the system 
is locally commensurate. It is our understanding that at 
long length scales the system is actually incommensurate 
with discommensurations (or kinks) separating the lo- 
cally commensurate regions. On length scales long com- 
pared to the distance between kinks, this state behaves 
like an effective "elastic solid" which in one-dimension 
has the same quantum critical behavior as a Luttinger 
liquid. Thus, this state is qualitatively equivalent to our 
weak coupling picture, albeit with substantially renor- 
malized parameters. 
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APPENDIX A: FERMIONIC HAMILTONIAN 



We considered an extended Hubbard model on a ladder with a Hamiltonian of the following form 

H = -tJ2 { C ha C i+hjo + h.C.} -tj_J2 i C i,0a C iM + h - c -} + U Y1 n i^ n i,ji + V \\Y^ n i,3 n i+hi 

(io) i,j i,3 

+ V±~y] n lti n 2 ,i + V d y^(ni,inj + i,2 + ^,2^+1,1) + J±^2 ■ S ii2 + J\\ S itj ■ S i+ ij (Al) 

i i i ij 

Here i labels the sites along legs and j labels the legs (or rungs); the coupling constants U, Vy, V± and V d represent 
the on-site Hubbard interaction and various nearest and next-nearest neighbor Coulomb interactions, and J± and Jy 
are the Heisenberg interaction along the rungs and chains respectively. 

After diagonalizing the kinetic part, we can rewrite the above Hamiltonian with the right and left movers of 
the bonding and anti-bonding bands represented by the operators ipRi,ipLi,ipR2,ipL2 as below, where ipi l2 (x) = 
(cji ± Cj 2 )/\/2a. In the low energy limit, the free part of the continuum Hamiltonian density can be written as 

7T 2 

^o = v fl {-(J L1 J L1 + J m J R1 ) + -ir(J L1 ■ J L1 + J R1 ■ J R1 )} 

7T 2 
+ V fl{ 2 ( J L2 J L2 + JR.2JR2) + 7^(J L2 ■ JL2 + J R2 ' ^2)}, (A2) 

where Ji.R.L — ip\ R ^u^i.R.Lcn and Ji.R.L = ip\ R L a ipi,R,La, are the right and left moving components of the charge 
and spin current density for the bonding (i = 1) and anti-bonding (i = 2) fermions. 

The interaction part of the Hamiltonian splits into several terms. First, we have a set of terms involving only the 
charge currents: 

'Hvnt.c = + -^{V\\ +Vd) + —^{JirJir + JilJil + J2RJ2R + J2LJ2L) 

+ {j + Wl + W - ^y^) + ^ - § cos2 fc/1 J„ - ±J ± }j 1R J 1L 

+ {j + (Vj| + V d )(l ^|^) + I cos2fc /2 ^, - A J ± } J 2K J 2i 

\J 1 1 3 3 3 "1 

+ {- + Vy(l - -cosfc_) + V d (l + -cosfc_) + -V±_ - -J|| cos/j_ + — Jj_|(JiflJ2fl + JilJil) 

+ { j + Vj|(l - i cosfc+) + V d (l + i cos/c+) + I VI - |j|| cos/j+ + Aj ± |(j 1K j 2L + J 1L J 2R ), (A3) 

where fc+ = fc/i + fc /2 = tt(1 - (5), = k fl - k f2 = 2 8m" 1 [t±/(2 cos ^)] . 
Next we have the couplings involving the spin currents: 

Wmt, s = { — — + -y- + -£-}(JirJir + JlLJlL + J2RJ2R + J2LJ2L) 

- {u + 2(V[| + KO cos 2^ + VI - I J±}JikJil - {t/ + 2(Vj, + V d ) cos2fc /2 + V± - | J±} J 2fl J 2L 

- {U + 2(Vj, - V d ) cosfc_ - VI - (1 + i cosfc_)J|| - ^}(J lfl J 2fl + J 1L J 2i ) 

- {[/ + 2(V], - Vc) cosfc+ - VI - (1 + 1 cosfc + ) J,, - ^}(J lfl J 2L + J 1L J 2il ). (A4) 
Next we have the low energy couplings associated with singlet-pair and triplet-pair tunneling: 

H int ,pt = {f/+(2(V|, - V d )- ^J||)cosfc /lC osfc /2 - V ± + ^J ± }{A\A 2 + h.c.) 

+ {2(V]| - V d ) + ^}sinfc / i sin k f2 {A\A 2 + h.c.) (A5) 

where A = (tp R ]-ipLi — V'-RJ.V'lt)/'^ ^ s the singlet pair operator on a given chain and A is its triplet counterpart. 
Notice that 1 and 2 stand here for the chain label. 
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Finally, the low energy Umklapp scattering terms are 



2 Hi 



{{j + j**[V\\(\ cosfc_) - V d {\ +cosfc_) + |j H ] + ^Vi + ^) N M 



- {U + e^{-2{V\\ -V d ) + ( l -+ cos k_)J\\)-V±-'^-}NiNl 

+ {^-(V\\-V d - \j\\)e i5 * - ^ + | Jx}(mt li? m 2i + m^mu)} + /i.c, (A6) 

Here = ip Ra ipLa and = ip RrT (a/2)ipi J(7 are CDW and SDW (Neel) order parameters, respectively, m is the 
paring order with 2k f momentum, for example, ttlr — ^Pr.^r.i- 

Following the standard Bosonization procedure with the assumption of Eq. 12.11 we arrive at the bosonized Hamil- 
tonian density in the section II. The bare values of the weak coupling constants are given as 

3 

g c+ = U + V\\ [4 + costtS (1 + cos £;/_)] + 2 V± + V d [4 - cos ttS (1 - cosfc/_)] + -,/ cos7r6(l + cosfc_), 
3 3 

,9c- = — (Vj| + - J|j) COS7T(5 (1 — COS kf-) — V± + COS7T<5 (1 + COS kf-)V d — — J_L, 

.g s + = U — Vj| cos7r5 (1 + cos kf-) + V d costtS (1 — cos kf-) — ^-(1 — — cos7r<5) — 

g s - = Vj| cosnS (1 — cosfc/_) + V± — V d cosnd (1 + cosfc^_) + -^(1 ~ — cosnS) — -Ji, 

g 3 = 2(Vf| -V d + -jt)[cosk- +costt5], 

1 3 

,g 4 = f7 + 2(Vf| - Vd)cosfc_ - Vj_ + J||(costt5- -cosfc_) + jJj_, 

5uc = U - 2costt(5 (V|, - Vrf - - VI + ~J X , 

5 „ 3 = t/-2cos7r5 [VJ| — — Jy(j + icosfc_)] -Vj.-^, 

,g U4 = C/-2cos7r5 [(V||+y d )cosfc_- J||(l + icosA;_)] +V1 + ^, 
.9i = 5s+ + 3s-, 32 = 3.s+ - 3s-, 35 = 34 - 33, 3«5 = 3« 4 ~ 3»3, (A7) 

where cosfc_ = 1 — t 2 L /(2t 2 ). Up to the first order, these coupling constants are independent on the doping 6. 

When away from the half-filling, the particle-hole symmetry Eq. 12.11 only holds approximately at small doping 8 
as k+ — 7r = 6 7T, Avf /a = S tj_ir. Taking these into account, there are some small residue terms as appearing in Eq. 
12.61 12.71 12.81 they vanish linearly with doping. The corresponding coupling constants are 

13 1 

A 3c = {77(^1 + Vd) + o ^||}sin7r<5 smk f _, Ag s = --(V« + V d ) sirnrS sinfe/_, 

X I COS k 

Ag uc = -2 mutt* (V], -V d - - J,,), Ag u3 = -2 811171(5 [V], -V d - J||(- + — g^)] 

Ag u4 = -2sin7r5 [(V|| + Vd) cosfc_ - J||(l + i cosfc_)] , Ag u5 = Ag u3 - Ag u4 . (A8) 



APPENDIX B: BOSONIC REPRESENTATION OF THE ORDER PARAMETERS 

The difference of the charge density between two legs reads J2 a {—y +1 Cj a (i)cj a (i) = J2a V'LWV^o- (%) + 
tpl a (x)ipicr{x)- After expressed by the right and left movers, it contains the staggered part, ie. Ocdw ■ Similar situation 
happens to its triplet counterpart OsDW,z,x,y Using the bosonization identities: tPr,l{x) = l/V2na exp{±iy / 7r(0(a;)± 
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0(x))} and we can obtain their bosonic expressions as below, 



Ocdw{x) 

OsDW,z,x,y(x) 



2Y 



i'lLo ( x )i>2Ra (x) + 1p\ La (x)ip 1Ra (x) + h.C. 

^\ La ( x )^/ 2 )^2Rfi(x) + i>\ La {x) {a /2) a pip 1Rf3 (x) + h.c. 



oc — { cos(v / 7r0c+ — Sirx) 
ira 1 



+ sin(v / 7T0 c + — Sirx) < 



2 cos \pKQ c - cos v / tt0s+ sin ^/nO s 

{COS ^/TT(j) s+ cos ^/n9 s - 
COS Y / 7T0 S+ COS ^/TT(j) s - 
— sin t/tt9 s+ cos ^/tt4> s - 

—2 sin ^/ir9 c - sin y / 7t0 s+ cos ^JHQ s ^ 
sin y / 7T(/) s+ sin y^s- 
— cos y/TT0c— ^ sin ^/n9 s+ sin y^s- 
cos \/t^0 s+ sin y / 7T(/) s _ 



(Bl) 



where T equals i?7j(l)?7j(2) for the singlet and ^-component of the triplet order parameters, and ir/j(l)?7j(2) for 
x, y-components of the triplet order parameters, and the same as below. 

The difference of the bond strength between two legs is ^2j a (—y +1 Cj a (i)cj a (i + 1) + h.c. = ip\ a .(x)ip2a(x + a) + 
tpl <J (x)ipi cr (x + a) + h.c, similar is its triplet analog. Their staggered parts, Osp and Og P 



Osp(x) 

OsP,z,x,y( X ) 



= (-) x 2sm(k fl + ?-6)i{e- iv6x - i5v / 2 



^l La ( x )(^/ 2 )a^2R0(x) + tpl La {x)(a/2) af3 ip 1Rfj (x) - h.c. 



2T _ 

oc — { cos(V7T0 c + — Sttx — Stt/2) < 
ira 



+ sin(v / 7T0c+ — Sttx — Sn/2) < 



2 sin yfitOc- sin yjit<j) s+ cos ^/tt6 s 
sin y / 7r</> s + sin ^/nd s - 
cos ^/tt9 c - < sin y/n6 s+ sin \fn^> s - 
cos ^/tt9 s+ sin \pK$ s - 

2 cos \fiid c - cos v / 7r</>s+ sin y/n9 s - 
COS y / 7T(/> s+ cos ^/ttQs- 
sin y/ndc- < cos y / 7r# s+ cos -y/^s- 
— sin v / 7T#s+ cos y/n<j) s - 



(B2) 



It is clear that Ocrw an d Ogp are real and imaginary part of i^ 1La tp2Ra + V4,LcrV'i-Rer respectively. 

Next we present the staggered part of the diagonal current density, i ) J ' +lc j W c j+i(* + 1) — ^- c -i an d its triplet 

analog as below 



^bcW 1 = (_)x sin ^ fc , 1+ a 7r /2){ e -*^-**/2» J ^LWteW -^ iff (a;)Viij CT (a;)) + /i.c. 

^ DC,z,x,y 



2r r / ^ x <K 

CX •{ COS(V7T'?>c+ — 7T<JX — )< 

7ra 2 



+ sin(v / 7T0 c + - tt<5x - -^) < 



2 cos \pxQ c - sin y/TT(j) s+ COS y/Tr8 s - 

!sin y / 7r0 s+ sin y^tf 
sin y / 7t0 s+ sin 
cos ^0 S+ sin y^ 

2 sin ^JttOc- cos y / 7r0 s + sin i/7r# s _ 

COS y / 7T(/) s + cos ^/ir6 s - 
cos s/kOc- \ cos \Pk9 s+ cos ^/lT(j) s - 
— sin \/tt9 s+ cos y/w(j) s - 



(B3) 



The difference of the current density along the legs is (c^ (?)c Jcr (z + 1) — /i.e.). Its staggered part is 



(57T 



(-)* 2cos(fc 1/ + -)(-i) ^ {e-^-^/^^^^^flW - VUWV'i^W) - /i.e.}. 
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Similarly, the staggered current along the rung i(t4«r(*) c io , (*) — h.c.) is: 



It can be shown that they satisfy the continuous relation |36j| . so does its triplet counterpart, Oddw* staggered 
currents along legs and rungs have the d-wave feature. We use currents along the rung as order parameters. Their 
bosonized forms are 



Oddw{x) 

ODDW t ,z,x,y{x) 



2r 



oc — { cos{\fTT 4> c + ~ it 8 x) 



+ sin(v / 7n/>c+ — ir5x) 



2 sin \ZttO c - cos y / 7T0 s+ sin ^/tt9 s 

{COS y/n(f> s+ COS y/lT0 s - 
COS V / Tr0 s+ COS y/TT<j) s - 
- sin ^/tt6 s+ cos y/n</> a - 

2 cos yfizQ c — sin i/7T<^ s + cos y^^s 
sin v / 7r0 s + sin y/Tr9 s - 
sin y/ndc- ^ sin ^/ttQ s+ sin \pK^ s - 
cos y/nOs-]- sin y/ir(f> s - 



> . 



(B4) 



It can also be seen that the DC and DDW order parameters are the real and imaginary parts of ip 1La ip2Ro- — '4 ! 2La' l l )l R'*-> 
respectively. 

Finally, the bosonized forms of the d-wave and s-wave pairing order parameters are 
2^(1)^(1) 



e iV¥0 c+ ^_ cog ^q c _ sm y^f0 5 . + s in y / 7t0 s _ + i sin \fiiQ c - cos yfK(j) s+ cos \fir(j) s J) 

A s = (^lif^liU _ V'lI^lflT) + {lp2L\i>2Ri ~ 1p2Li^2RT;) 

_ V]{ )Vl{ ) e i^Re c+ ^ cog ^q c _ cos y^f0 s+ cos y/T:(j) s _ + i sin \/t:6 c - sin v / tt0s+ sin v / tt'/'s-) 



(B5) 
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